No production of entropy in the Euler fluid 
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Abstract 



OO \ We derive the Euler equations as the hydrodynamic limit of a stochastic of a hard- 

sphere gas. We show that the system does not produce entropy. 

o 

; 1 Hydrostatics of a gas of hard spheres 

We take space to be A C (aZ) 3 , and suppose the length a, representing the diameter of a 
molecule, to be so small compared with the variation of the macroscopic fields that we can 
replace all sums over A by integrals. The possible configurations of the fluid are the points 
in the product sample space 

cceA 

so a configuration is specified by the collection {cox}xeA- For each x, 



P : L J 

Here, e is a small parameter having the dimension of momentum. If the system is in a 
configuration oj, such that ujx = 0, then we say that the site x is empty. If u>x = k, we say 
that the site x is occupied, by a particle of momentum k. This simple exclusion of more 
than one particle on each site incorporates the hard-core repulsion between the particles, 
^ . which are thus hard spheres sitting at some of the points of A. 

The state of the system is a probability on ft, denoted by \i. We denote the set of states 



by S. The 'slow variables' of our model are the 5 extensive conserved random fields 



if u x = 

1 if lo x = k 



£x(oj) 



[OJ] 



if ujx 

k • k/2m + $(x) ifu! X 

if uj x = 
k if uj x = k 



Here, $(a;) is the external potential energy per particle. The slow variables appearing in 
hydrodynamics are the n = 5|A| means in the state \i: 

N x = E At [Mx\, E x = [£ x ] \ ^ x = [V x ] . 
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The (von Neumann) entropy of any state /x is 

S(p) :=-k B ^2n(w)]ogn(u). (1) 

In information geometry, the choice of slow variables {X\,X2 ■ ■ ■ ,X n } defines the in- 
formation manifold A4, which consists of states of maximum entropy among all states with 
given means, say 

E^Xj] :=fi-Xj = rjj, j = l,...n; 

By the Gibbs-Jaynes principle, such a state has the form 

H(u) = exp | - €j x j H - !og H 

where the dual, or canonical, variables £j are Lagrange multipliers determined uniquely by 
the means r/j. In our case, for each x, the duals to energy, number, and momentum are, 
respectively, Px,£,x,Cx, an d so the state has the form 

H(u) = Y[E x 1 exp{-£xMx{u) - p x £x{uj) - Cx ■ V x {uo)} . (2) 
x 

Such a state is said to be in local thermodynamic equilibrium, LTE. Equilibrium holds when 
(3... are independent of x. In [7] we found the (grand) partition function at point x, 

( 2itm \ 3 / 2 

Ex = 1+ Wfa) 6XP ^ x ~ + m C» ■ Cx/Wx} 

by replacing the sum over the momentum lattice of size e by a Gaussian integral. The 
product structure of an LTE state means that an observable at a point of A is independent 
of an observable at any other. Note that Hx{@} = (1 — Nx)- It then follows from (J2J) that 

l-Nx = S~ 1 . (3) 

If w x 7^ 0, the state // can be written in Maxwell form 

//(sb, k) = N xP (x, k) = NxZ- 1 exp{-px$(x) - f3 x k ■ k/(2m) - Cx • k} 

where Z is the canonical partition function: 

I?* J + -55-}- (4) 

We note the identity for each x 

where Zq is the canonical partition function when $ = 0. The external potential does not 
influence the local velocity distribution p, as it is cancelled out by the partition function. 
The mean fields are related to the canonical fields by 

d i - a TV \ ( \ , 3 m Cx • Cx 



li « - -"^riogr.a; 



50 ° * A 



a; 
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The hydrodynamic variables are the mass-density p = mNa" 3 , the velocity field u = 
— C//3, the temperature = (k B /3)~ , and the potential energy per unit mass, <f> = $/m. 
We shall therefore eliminate £, f3, C, 3 in favour of p, 0, v? . Only £ remains to be found: 

e-te =Z X 1 N X /(1-N X ), (5) 

or 

£ = -p$ - log iV + log(l - AT) + log Z (6) 
where Zo is the canonical partition function when = 0. We need its gradient: 

„ _ 35,9 m uV ^ m«* ,■ d<N djN &•$ $ 9,6 

^ = ^e--2A-^ + ^ "iv--r^v-Aie + A-^- (7) 

We denote by iV and II the total values of the mean energy, number and momentum; 
then Q gives for the entropy at equilibrium 

QS(jt) = E + k B 0£N-u • n + A; s GlogS. (8) 

Compare this with the thermostatic formula 

OS = E + k B Q£N -u-U + PV (9) 

where P is the pressure and V is the volume; we see that 

P = k B ^eiog~ = -k B ea~ 3 log(l-N). (10) 

2 The fundamental equation 

Our model for the dynamics is a path {p-(t)}t>o through S(f2) determined by the Liouville 
motion, interrupted by a thermalization at random points x, occurring at random times. Af- 
ter a thermalization at x, the state changes to one in which observables at x are statistically 
independent of those at any y 7^ x, and, restricted to Six, the state is in thermodynamic 
equilibrium. The law for the random time between collisions is of exponential form, but 
the rate of the process depends on the local density of the gas, and its temperature, and so 
on the state itself. This means that the dynamics falls outside the usual theory of Markov 
chains, in which the updated state is linear in the current state, and the transition matrix 
is the same for all states; we are in the business of non-linear Markov chains. Nonlinearity 
itself is not the problem; the ultimate goal of this work is to put an external potential $ 
into our version [7] of the Navier-Stokes equations, which are non-linear. However, we are 
going to use the concept of conditional probability to derive the master equation, and care 
is needed if we stray from a path in Vl to a path in S(O) which does not come from a path 
in $7. We shall adopt the following way out, which makes use of the assumption that the 
density is low, and so the collision probability is small. In calculating the probability per 
unit time that a configuration uj at time t S (to — t, to) move to another point of O, we 
shall assume that the hydrodynamic parameters satisfies the Euler equations, (|3T)1) - (|3*21) in 
the small time interval (to — T,to)- The initial values of the hydrodynamic parameters in 
these equations are taken to be those of the true state at time to — T. Here we only need to 
consider T < 20ti, since the survival probability falls exponentially. We then assume that 
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the Markov process between to — T and to is linear as usual, and takes place in an ambi- 
ent Euler fluid which determines the rate of thermalization. In this paper, we show that, 
neglecting collisions, the means of the slow variables do indeed satisfy the Euler equations, 
showing self-consistency. We can then show that there is no production of entropy in this 
case. 

The Liouville dynamics is determined as follows. If at time t = a configuration uj £ 0, 
has a particle at x, then it follows Newton's laws 

dx k dk „ . , 

— = -, =_V$:=m/ 11 

at m dt 

for a time. Before we take the continuum limit, space is A, which is discrete; a particle 
following Hll[) will nearly always leave the lattice. In that case, after any given random 
time t after which a thermalization occurs, we place the resulting thermalized particle at 
the lattice site nearest to x(t), ties being decided by tossing a coin. In the limit o->0we 
expect this to introduce a negligible error in the location being assigned to the particle. 

Let (i £ S(O) and define Nx = ^{^x 7^ 0}> the probability that x is occupied in the 
state \x. Then N x = [J> ■ Nx- Let px(k) be the conditional probability that ujx = k, given 
that N x = 1: 

Px(k) = n{ui : ui x = k\ui x ^ 0} = N^fifa : ui x = k}. (12) 



Thus, 



n{ui :ujx=k} = N x p(x, k). (13) 



This does not mean of course that Nx and "P x are independent in the state \i. 

Given that no collisions occur, the particles obey (|11|). and this induces the Liouville 
motion on the states, namely, Boltzmann's equation with no collision term: 

M^L + A . a (a . 5 k) + fm Vkfl{x ^ k t) = . (14) 

ot m 

This dynamics of the probabilities does not correspond to an underlying motion in O for all 
time. We can find initial conditions for which two particles both arrive within a distance 
a/2 from the same point at the same time. Such a configuration does not lie in O. Our 
assumptions of no collisions is true with high probability for a few time steps, but colli- 
sions are almost sure to occur eventually. We will replace the problem of tracking which 
collisions actually occur if we follow (|11[) by Boltzmann's idea of introducing a probability 
for collisions. Instead of giving, as Boltzmann did, the probability that a pair of particles 
with given position and momentum produce another specified pair, it will be enough to 
assume that the collision is 100% efficient (in the terminology of [0]). This means that after 
a collision at x, of a particle with momentum k, we replace this configuration by a particle 
at x with momentum k', with probability determined by the Maxwell distribution of mean 
momentum k and mean energy k • k/{2m) + By construction, this does not alter the 
means of the slow variables. We shall refer to this event as a thermalization, rather than 
as a collision. We shall assume that the mean free path is large on the molecular scale, 
and neglect the possibility that a snapshot of the lattice catch a particle in the process of 
thermalizing: with high probability it will be in Newtonian motion between collisions. From 
any initial configuration uj, we can follow the dynamics as time progresses, following this 
dynamics. It is a random path through O, that is, a process, defined for to — T < s < to, 
and described by the family of points {/j,(s)}t -T<s<t - 
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Let t denote the random time between collisions, whose law is determined by the condi- 
tional probability w(x, k,to;t)dt that a particle, certainly at x at time to with momentum 
k, travel under Newton's laws a free time t and thermalize in the interval (to +t, to + t + dt). 
Since it will thermalize eventually, 

/•OO 

/ w(x, k,t ;t)dt = 1. (15) 
Jo 

The important small parameter is the mean free time, also called the relaxation time, tf 

ti(x,k,to) ■= / tw(x,k,to;t) dt. (16) 
Jo 

This is invariant under Q, unlike the mean free path, which is related by £ = kt^/m. For 
free particles, = 0, these were found in [7J. They are not affected by the potential, so we 
use the same values here. We also introduced 

C(y, k,t)dt = Probja particle at (y, k) at time t thermalize in (t,t + dt)} (17) 
W(x, k, to; t) = Prob {a particle at (x, k, to) survive up to time to + t} . (18) 

Put z(t) := (x(t), k(t)) := Tt(z), the solution to Newton's equations with initial point 
z := (x,k). Then we have the relation 

w(x, k, t ; t) = W (x, k, t ;t) C (z(t),t + t ) (19) 

and just as in [7] we get from (jTo^l . 

W(x,k,t ;t) =exp|-^ C{z(t x ), t + h) dtA . (20) 

As in |7j, eq. (36), C is proportional to the density p. Thus if p is bounded below by a 
positive number, we get an exponential decrease for W, the survival probability, along an 
orbit. 

The fundamental equation will relate the probability distribution p,(z, to) at an arbitrary 
point z in phase space at the time to to a Maxwell distribution /I at the same point, with 
small correction terms. Let z(s) denote (for s > 0) the point in phase space on the backward 
Newtonian orbit of the point z at time to — s. For each sample path uj(-), in which z is 
occupied at time to, there is a unique time to — t at which the last thermalization occurred, 
in that no further thermalizations took place on the orbit {z(s)}o<s<t- It follows that the 
free time of the particle thermalized at to — t must have a free time t! say, with t' > t. 

In the following, 'event' will mean an event in the sample space (the path space of f2) 
of the process constructed above in the time interval (to — T, to). Let E(s) denote the 
event, that there is a particle at z(s) at time to — s. We shall consider the conditional 
probabilities of collision and free motion along the orbit, given E(0). We shall then recover 
a formula for p(z(0)) using Bayes's theorem. Liouville's theorem, that d s xd s k is invariant 
under Newton's laws, enables us to deduce equations relating the density of probability 
from equations relating probabilities. For example we shall write Pr{£'(0)} = p(x, k,to), 
referring to the densities. 

Let F(t) be the event, that ^(0) occurred, and the last collision occurred at time to — t 
producing a particle at the phase point z(t). Then F(t) C E(0), since the event produced 
exactly the right state, z(t), which evolves to our point z(0), as it undergoes no further 



5 



collisions. Let F(a, b) denote the event, that the last collision was at some s, a < s < b. If 
(a, b) n (c, d) = 0, then F(a, b) and F(c, d) are disjoint. So, assuming smoothness, F(t) has 
a probability density on R + , say f{t). Let G(t, t') be the event that z{t) is occupied at time 
to — t, and that the particle's free time is t'. Let H(ti) be the event that F(s) occurred for 
some s < t±. We shall need the formula 

Pr{F(t)}/Pr{tf(t')} = /(*)/ f f(s)ds 

Jo 

which is true if / is smooth and /(0) is not zero. We see that G(t, t')PiH(t') C E'(O). Hence 
certainly 

G{t,t')r\F(t)r\H(t') c E(0). 

Hence 

Pr{G(t, t') n F(t) n tf(t') n E(o)} = Pr{G(t, t') n F{t) n tf(t')}- 

Thus 

Pr [G{t, t') n F(t) n H(t')\ E(0)} := Pr{ J E(0)}" 1 Pr{G(t, t') n F(t) n F(t') n £(0)} 

= M (^(0), to)" 1 Pr{G(t, t') n F(t) n F(t')}. 

Let Pr#/{.} denote the conditional probability of an event, given H{t'). Then from what 
we have just seen, 

Pi{G(t,t')nF(t)\E(0)} := Pi{G(t,t')nF(t)nH(t')\E(0)}PT{H(t')}- 1 

= Pr{G(t, t') n F{t) n Hit')}^- 1 Pi{H(t')}-\ 

Now, 

' dt' J dtPr{F(t)nG(t,t')\E(0)} = 1. 
as the last collision must have happened for some t and some t' > t. Hence 

roo rt' 

/ dt' dtPr{G(t,t') nF(t) nHit'^^PiiHit')}- 1 = 1. 
Jo Jo 

As n(z(0),to) does not depend on t or t', we get 

uum t \ - f°°dt' ^{G(t,t')nF(t)nH(t')} 

^(o),t )-y o dtj o dt • 

We can omit H from the numerator, as t < t' is enforced by the region of integration. So 

= r dt , r'/^on^) 

^ Jo Jo Pv{H(t')} 

f°° df , f , Pr{G(t,t')\F(t)}Pr{F(t)} 
~ Jo Jo Pr{H(t>)} ■ [22) 

Also, if t < t', then 

Pr {G(t, t')\F(t)} = Np(z(t),t - t)w(z(t),t - t; t'), (23) 
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(what we have been calling a collision is a thermalization.) In applying (|21j) . we note that 
the case /(0) = corresponds to no collisions, so we may assume that /(0) 7^ 0. We expand 
(|2M|) up to first order in i, as i < t', and t' <T = 20tg = Oitg); this gives 



Pr {G(t, t')\F(t)} = Jl(z, t )w(z, i ; t') 
Put this, and use (JUJ), in (J22|) to get 







ji(z,t )w(z,t ;t') 



dt' 



H(z,t )w(z,t ,t') 



m 



(24) 



The unknown term involving /'// does not appear, because Jq (t — t' /2)dt vanishes, and 
Jq (t — t'/2)tdt is of second order in t' , and so can be neglected. We can now do the integral 
over tf, using ()15j) . and ((TB|) . This gives 



H = fi 



m 



_d_ 



(25) 



This equation was derived in |7j using another method (aka guesswork), in the case when 
$ = 0. 

It does not seem fruitful to seek accuracy up to and including 0(tf). This would involve 
introducing unknown parameters, such as /'//; worse, we would have to keep terms of 
second order in taking the continuum limit of the lattice; this introduces diffusion terms 
into the equations, which come from the Ito corrections to calculus. A similar extension 
of the work of Chapman and Cowling [2], who start with the Boltzmann equation rather 
than a master equation, is generally agreed not to have been worth the effort. Keeping 
only terms of first order in tg leads to the surprising result that the equations of motion 
are determined, knowing only that fx is in LTE. We do not need to assume, as is done in 
PQ, that the means of the slow variables in the approximating LTE-st&te ~p are the same 
as the true means, in the state fi. Indeed, this turns out not to be the case; the reason 
is that Ji is conditioned by the fact that a thermalisation has occurred, and is not just 
the maximum entropy estimate of \x. Then (J,(x, to) is the sum of such terms from various 
nearby points (x — kt/m,to — t) and for a given fc, all the contributions are from one side 
of x, so the means should differ unless the state is constant in space and time. In this 
paper we neglect the difference between \x and /J, which is the cause of dissipation in the 
Navier-Stokes equations. We arrive at the Euler equations, and show that they are free of 
dissipation, as expected. 



3 The Euler equations 

The Euler equations follow from the approximation of zero th order, in which the difference 
between ll and JI is neglected. Mean fields for states in LTE are computable in terms of 
Gaussian integrals. 
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The velocity field of a particle is the random field Tx '■= T > x/ m , arid the mean current 
of a random variable %■> conserved or not, is, in the continuum, low-density limit 



J x = J d 3 kpT X . (26) 

Our assumption is that on thermalization there is no change in the means of J\fx, T'x or 
Ex- Since the space integrals of Nx and Ex are to be conserved in time, their equations of 
motion are of the form 

P + d d 4 = (27) 
E + dj4=0. (28) 

Momentum is not conserved; in a small volume at x a particle of momentum k in time dt 
is changed to k — V<3? dt, so the momentum obeys 

^ + a.jJ, + ^ = 0. (29) 
dt 3 m m 3 v 1 

In [7] we showed that we can evaluate the currents ()26|) of hydrodynamics exactly if p is in 
LTE, with parameters p, u l , say. By the same method, (|2*7|) - ()29|) can be written out: 

^ + d j Qm J ) = (30) 

(pui) + dj (puiUj) + diP + pdi4> = (31) 
f3k B „ , 1 , a f /3fc 



1/ V 2m 2 J J m 



9, {/m, Q^-G + P/p + ^ W + A J = 0. (32) 



Here, = <I>/m. Note that the pressure P appearing in these equations is that of an ideal 
gas, not the equilibrium pressure of the model; the two differ by terms of second order in 
the density. This reflects that the derivation of the fundamental equation depends on the 
assumption that the density is small. 



4 Entropy production 

From von Neumann's formula, 

S:=-5>Mlog0(") (33) 

UJ 

we conclude that 

s = - £ AM log //M - £ AH- (34) 

UJ UJ 

The last term is zero, so we may regard — log p as a random variable S whose average rate 
of change is given in the Schrddinger picture by p ■ S. The current of the entropy would 
then be, per site 

a 3 J>(x) = S(u)T j = £NT j + (3ST j + CViT j + log 3 T j , (35) 
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when [i = Np and p is in LTE. Its mean density, J|, in the state p has three terms: 



(a~ 3 p ' f A/"T J ) + a~ 3 log Eu • T = - log(l - iVW 

a- 3 CV-(P 4 T^) = C^4 



a~ 3 ' 



We expect that entropy will be conserved, because of Kossakowski's argument [3]: in the 
absence of collisions, the projections onto the information manifold do not create any en- 
tropy if time is continuous. See also 0. Without collisions, Boltzmann's equation also give 
zero entropy production, if the assumption of LTE is made. So we expect that the rate of 
change of entropy is balanced by the outflow though the boundary: 

s + 3j + CJi< + PJ> B - log(l - N)u^ = 0. (36) 

where s denotes the entropy density. On the other hand, we have 

a 3 s = ZN + CILi + f3E + log 5, 

giving 

a 3 — = Nd^ + Ihd t C + Ed t p + -(^d^ + ^d t C+^d t (3) (37) 



dt E\d£ dQ 8(5 

HdtN + Cd t Ik + (3d t E. (38) 

Now, 

1 5S 

-7^7 = -iVetc. 
so the first line (|37|) vanishes, to leave 

a 3 d t s = &N + Cd t Ui + 08 t E. 
We divide by a 3 to get the densities p, a~ 3 E and Wi, and the relation 

/ 8 t sd 3 x = [ (^dtp + CdtVJi + (3dtE/a 3 )d 3 x (39) 
J A JA\m J 

= ^(Cpfi + drfji+dtfJit+drfJi) (40) 

- | A (^ + (te f + WW+^)^ (41) 

where n denotes the normal direction to the boundary. The term in (|41|) involving /j = 
—di<fr/m represents the free-energy change due to the work done by the external field on 
the particles; we shall see that this is cancelled out by another term, showing that this work 
done is not thermalised by the dynamics. In spite of the apparent Onsager form of (|40[). 
we cannot identify it as the entropy production and (|41|) as the flow through the boundary, 
because the surface term differs from the entropy current, as in (|36|) by — log(l — N)u 3 . 
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Indeed, if we put the Euler currents in (|40[) we do not get zero, even when = 0. We must 
therefore add and subtract the missing term, to get 

8 t s = J (Cpf t + d&Ji + i)jCH, + djpJi - djia" 3 log(l - N) Uj )) d 3 x (42) 
Jidai. (43) 

9A 

If we now put in the Euler currents, (|30j) - (|32() . we still do not find exactly zero for the entropy 
production, (|42j) . Indeed, the entropy production involves the logarithm, whereas the Euler 
equations are polynomial, so cancellation is not possible. We find that (|4*U|) vanishes up to 
terms 0(N 2 ), that is, in the limit of low density. This reflects the low-density assumption 
inherent in the axiom that the BBGKY currents, l|26[). are the actual currents of particles 
flowing out of the region d 3 x. While this is a reasonable model for point-like particles, it 
neglects the fact that if the lattice-site just beyond the boundary is occupied, then a particle 
moving out of the region will not be able to land, as no configuration in the sample space 
can have two particles at the same site. Our dynamics did not need to specify the rule as 
to what will happen, as (in the low density limit) the probability that the point is occupied 
is small. 

It turns out that if we modify the pressure in the Euler equations, to be the equilibrium 
pressure of the gas with hard core, (|10|) . rather than that of the ideal gas, then the entropy 
production is exactly zero, in conformity with Kossakowski's theorem: 

Theorem 1 Take the pressure P in 1)31)1 and (|32|) to be that of the hard-core gas at 
equilibrium, (jTU|) . Then the entropy production, is zero. 

PROOF. We first recall the relation between the canonical variables and the hydrodynamic 
variables. We saw that 

1 iii 

P = 1 - R C 3 = -f3u 3 



k B e * k B e 



so from (JJJ) we get for the entropy production (|42|): 



•Sprod — C Pf'i 



pu J / 3djQ m u l u l m u % i djN djN m dj<p m 4> a n 

+ (P^ + S tJ P)(-^ + -^\ 



2 



U. 



+ a^Y 1 ^^ ~ a ~ 3 log ^ ~ N ) d i uj 
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